METHOD FOR REMOVING PSEUDO SIGNAL AND PROGRAM FOR 
REMOVING PSEUDO SIGNAL 



CROSS REFERENCE TO RELATED 
5 APPLICATIONS 

This application makes reference to, incorporates the same herein, and claims all 
benefits accruing under 35 U.S.C. § 119 from an application for Method For Removing 
Pseudo Signal And Program For Removing Pseudo Signal earlier filed in the Japan Patent 
10 Office on September 2, 2002 and there duly assigned Ser. No. 2002-256871. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

15 The present invention relates to methods for processing signals, and more particularly to 
a method for removing a pseudo signal having a lot of energy contained in the signals. 

2. Description of Prior Art 

A conventional method for processing signals using Fourier transformation is superior 
for grasping and processing a signal in the frequency domain, and are the techniques used 

20 in various fields. However, Fourier transformation is not suitable for grasping and 
processing a signal in the time domain in comparison with a method for processing signals 
by using wavelet transformation which has been researched and developed in several fields 
for these years. In the impact tests for such as satellites, when the acceleration signal 
measured contains a pseudo signal, since the prior arts have no effective way to remove the 

25 pseudo signal, the test data having frequent pseudo signals had to be given up for out of 
evaluation. 

Since the conventional method for processing signals using Fourier transformation 
presumes that the signals have periodicity, and since it processes the signals in the 
frequency domain, it cannot process the signal component useful as test data (hereinafter 
30 called "main signal component") by separating a pseudo signal components therefrom. On 
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the other hand, simple use of a prior art method for processing signals by wavelet 
transformation has a problem that even the main signal component can be removed when 
the pseudo signal component is not significant in comparison with the main signal 
component. 

5 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the present invention to solve the aforementioned 
problems, and the present invention provides a method for removing a pseudo signal from 

10 an original signal which includes the pseudo signal consisting of zero-shift components. 

The aforementioned problems can be solved by the present invention having the 
following features. The present invention includes a method for removing a pseudo signal 
from an original signal which includes a main signal free from a zero-shift component and 
the pseudo signal consisting of zero-shift components, the method comprising the steps of 

15 means for integrating said original signal with respect to time; means for decomposing said 
integrated original signal is decomposed into the sum of one Approximation (A(t)) and at 
least one Detail (D(t)) by wavelet transformation with a level number deepened until which 
level number reaches to a threshold level number where the ratio of the energy of said one 
Approximation to the energy of said original signal has a predetermined relationship to a 

20 predetermined threshold value; means for setting said Approximation to zero at said 
threshold level number as a pseudo signal component by a threshold operation; means for 
deriving an integrated main signal by reconstituting said Details up to said threshold level 
and said Approximation previously set to zero by wavelet inverse transformation; and 
means for deriving said main signal by differentiating said integrated main signal with 

25 respect to time. 

The present invention also includes that, in the method for removing a pseudo signal just 
described, said threshold level number is the minimum of the level numbers where the ratio 
of the energy of said one Approximation to the energy of said original signal becomes 
smaller than a predetermined threshold value. 
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The present invention further includes that, in the method for removing a pseudo signal, 
said original signal is an acceleration signal detected by a sensor for an impact test, and said 
pseudo signal is a zero-shift component induced to said sensor by the impact of said impact 
test. 

5 In addition, the present invention includes an apparatus for removing a pseudo signal 
from an original signal which includes a main signal free from a zero-shift component and 
the pseudo signal consisting of zero-shift components, the apparatus comprising means for 
integrating said original signal with respect to time; means for decomposing said integrated 
original signal is decomposed into the sum of one Approximation (A(t)) and at least one 

10 Detail (D(t)) by wavelet transformation with a level number deepened until which level 
number reaches to a threshold level number where the ratio of the energy of said one 
Approximation to the energy of said original signal has predetermined relationship to a 
predetermined threshold value; means for setting said Approximation to zero at said 
threshold level number as a pseudo signal component by a threshold operation; means for 

15 deriving an integrated main signal by reconstituting said Details up to said threshold level 
and said Approximation previously set to zero by wavelet inverse transformation; and 
means for deriving said main signal by differentiating said integrated main signal with 
respect to time. 

The present invention also includes that, in the apparatus for removing a pseudo signal 
20 just described, said threshold level number is the minimum of the level numbers where the 
ratio of the energy of said one Approximation to the energy of said original signal becomes 
smaller than a predetermined threshold value. 

The present invention further includes that, in the apparatus for removing a pseudo 
signal, said original signal is an acceleration signal detected by a sensor for an impact test, 
25 and said pseudo signal is a zero-shift component induced to said sensor by the impact of 
said impact test. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1 is a flowchart showing operation of the method for removing a pseudo signal 
according to the present invention; 

FIG. 2A is a schematic diagram showing the test system based on the method for 
measuring dynamic response characteristics by an impact acceleration meter used for a 
5 working example of the present invention; 

FIG. 2B is a schematic diagram showing the test system based on the method for 
measuring dynamic response characteristics by an impact acceleration meter disclosed in 
Japanese Patent publication for opposition purpose (Kokoku) No. 6-52270 used as a 
comparative example of the present invention; 
10 FIG. 3 is a diagram showing impact response acceleration signal a(t) having a pseudo 
signal; 

FIG. 4 is a diagram showing impact response velocity signal v(t) having a pseudo signal; 
FIG. 5 is a diagram showing main signal component D(t) decomposed by wavelet; 
FIG. 6 is a diagram showing a pseudo signal component A(t) decomposed by wavelet; 
1 5 FIG. 7 is a diagram showing the waveform of impact response acceleration signal a c (t) 

from which a pseudo signal component has been removed; 

FIG. 8 is a diagram showing the waveform of impact response acceleration signal a^t) 
of the pseudo signal component; 

FIG. 9 is a diagram showing the results of the evaluation for removal of a pseudo signal 
20 by Shock Response Spectrum (SRS); 

FIG. 10 is diagrams showing acceleration signal having a zero-shift and its SRS; 
FIG. 1 1 is diagrams showing acceleration signal from which a pseudo signal has been 
removed according to the prior art method and its SRS; and 

FIG. 12 is diagrams showing acceleration signal from which a pseudo signal has been 
25 removed according to the method of the present invention and its SRS. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

First, general notion of wavelet will be explained. Wavelet decomposition is a 
30 technique where a signal is decomposed by scaling (enlargement) and shifting (transition) 
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of a function called "mother wavelet" into the sequence of wavelet coefficients, each of 
which is a partial approximation, and the signal is expressed thereby. Such mother wavelet 
is usually a localized waveform (support compact). Therefore, the waveform of a signal 
having a violent fluctuation can be expressed in small scale, and the waveform of a signal 
5 having a gentle fluctuation can be expressed in large scale. Since a pseudo signal 
component such as a zero-shift generally tends to deviate to a signal having a gentle 
fluctuation, the zero-shift component can be distinguished and removed by using such 
tendency. A typical procedure for the process of wavelet analysis is explained hereinafter. 
First, a wavelet coefficient C(j,k) represented by the following equation is derived by the 
10 wavelet transformation of signal s(t). s(t) is generally represented as a sequence of time- 
series waveform sample value: 



15 where discretization scale is a = 2 j (j = 1, 2, 3, L (L is the level (hierarchical number) of 
wavelet)), shift is b = k2 j (k = 1, 2, 3, L), and ¥ is mother wavelet. In the present 
invention, by using Daubechies wavelet, the signal is separated by hierarchy according to 
wavelet coefficient C(j,k) as shown in the following equations: 





A,(t) 



D,(t) Level : J- 1 s(t) =A(t)+R (t) 



7"X 



A 2 (t) 



D2(t) 



Level :J=2 s(t) = A(t)+ £ D, (t) = A (f)+IX (r) 



A 3 (t) 



D 3 (t). 



Le.veJ..:J=3 s(r) = ^(0+ gA(0 



Level . :J=L s(t) = A(t)+ £ A (0 
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where 

k 

a = 2 J ' < 2 J 

5 Aj(t) represents signal components having gentle fluctuation and large scale. After 
calculating the energy ratio: 

10 between signal Aj(t) and original signal s(t), the wavelet transformation is executed with 
deepening the level number until the level number reaches to the level number (hereinafter 
called "threshold level number") at which the ratio of the energy E ; has predetermined 
relationship with a predetermined threshold value. Namely, energy ratio E ; is used as an 
evaluation standard for evaluating the level number required. Preferably, such threshold 

15 level number is the minimum of the level numbers where the ratio of the energy Ej 
becomes smaller than a predetermined threshold value. As a result, the minimum of the 
level numbers where the ratio of the energy E ; becomes smaller than a predetermined 
threshold value is determined as level J. The threshold value can be, for example, 99.7%. 
This is because the energy of A(t) having gentle fluctuation dominates most of the part of 

20 the energy of original signal s(t). Hereinafter this level J is called threshold level number. 
Aj can be regarded as a pseudo signal component. 

Next, a series of procedure in the specific numerical calculation method for the signal 
separation will be explained. 
1) Selecting type of wavelet 
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In the present invention, Db3 (Daubechies wavelet, order N = 3) was chosen. The 
following table shows the values of discrete mother wavelet function points which 
correspond to various values of N. 

5 Table 1 



Type 


Discrete mother wavelet function points (2*N) 


Dbl 


0.5 


0.5 














Db2 


0.3415 


0.5915 


0.1585 


-0.0915 










Db3 


0.2352 


0.5706 


0.3252 


-0.0955 


-0.0604 


0.0249 






Db4 


0.1629 


0.5055 


0.4461 


-0.0198 


-0.1323 


0.0218 


0.0233 


-0.0075 


Db5 


0.1132 


0.4270 


0.5122 


0.0979 


-0.1713 


-0.0228 


0.0549 


-0.0044 


-0.0089 


0.0024 














Db6 


0.0789 


0.3498 


0.5311 


0.2229 


-0.1600 


-0.0918 


0.0689 


0.0195 


-0.0223 


0.0004 


0.0034 


-0.0008 










Db7 


0.0550 


0.2804 


0.5156 


0.3322 


-0.1018 


-0.1584 


0.0504 


0.0570 


-0.0269 


-0.0117 


0.0089 


0.0003 


-0.0013 


0.0003 







2) Calculating low-pass filter and high-pass filter for decomposition and reconstitution 
The filters are calculated by using discrete points w(k), k = 1, 2, 6 of Db3. 

2N 

10 Reconstitution : Lo__R(k) = 0.707w(£) Hi __R(k) = (-\) k Lo _R(2N + \-k) 
Decomposition : Lo _D(k) = Lo _R(2N - k), Hi _D(k) = Hi _R(2N -k) 

3) Binary down-sampling 

As shown in the following equations, binary down-sampling is carried out by calculating 
15 convolution for signal s(t) and the filters: 
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cA x (k) = [Lo _ D(k) * s(k)] ■ sam _ - - • f cd, +1 (*) = [Lo _ * cA j (k)] • jaw _ </(*) 
cD, (*) = [/tf _ D(k) * s(k)] • sam _ - - - ( c£> /+1 (A) - [#/ _ * cd, (*)] - ram __ d{k) 

The binary down-sampling is carried out by using the following binary down-sampling 
functions: 

5 

x(k) = • sam _d(k) = y(2k - 1) , when the number of y is odd 
x(k) = y(k) • sam_d(k) = y{2k) 7 when the number of y is even 

This operation also diminishes the number of data into half. The following system flow 
10 chart shows the operation from level j to lower level j+1 : 



cA t 



Lo D 



Lo JXk)*cA(k ) F(k)sa m_d(k) 

—> cAj+i 



Hi D 



G(k) 



*2 



-> cD 



f+\ 



Level j HiJXk)*cA(k) G(k) sam_d(k) Level j+1 



4) Deciding level J of the constituting signal 
15 Level J of the constituting signal is decided according to the energy ration: 




For this purpose, the wavelet transformation is executed with deepening the level number 
20 until the threshold level number at which the ratio of the energy Ej has predetermined 
relationship with the predetermined threshold value. Namely, energy ratio E ; is used as an 
evaluation standard for evaluating the level number required. Preferably, such threshold 
level number is the minimum of the level numbers where the ratio of the energy E ; 
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becomes smaller than a predetermined threshold value. This is because Approximation 
representing most of the pseudo signal component can be obtained by the wavelet 
decomposition using such threshold level number. As a result, the minimum of the level 
numbers where the ratio of the energy E ; becomes smaller than a predetermined threshold 
value is determined as level J. In the present invention, since the pseudo signal component 
dominates the integrated signal energy, the threshold value for the energy ratio can be, for 
example, 99.7%. In this way, the pseudo signal component is distinguished. Namely, by 
removing Approximation of level J, the pseudo signal component is removed. 

5) Up-sampling 

The number of data is doubled by up-sampling which inserts zero between the data 
(between even numbers). 

6) Reconstitution of the signals 

Two signals are reconstituted by using convolution with reconstituting low-pass filter 
and high- pass filter. 

7) Constitution of the signal 

Signal is constituted by adding two signals formed in just described 6) together. The 
following system flow chart shows the operation from level j+1 to upper level j: 



cA 1 - 
Level j 



L/C 


Lo_R 


) M i 


+1 w 
f2 
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Hi _R 


f2 






< 



cA 



7+1 



Hi_rXk)*G(k) cD )+l (k) sam_u(k) Level j +1 
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These are the general procedure for wavelet analysis. 

Next, the method for removing the pseudo signal to which the above-described 
procedure is specifically applied will be explained. FIG. 1 is a flowchart showing the 
method for removing a pseudo signal according to the present invention. This method 
takes advantage of the characteristic that the pseudo signal component is emphasized by 
integrating the acceleration signal, and thus the energy ratio of the pseudo signal 
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component to the original signal becomes large. The method of the present invention can 
be carried out by using such as a computer. Namely, it is carried out by inputting the 
original signal to the computer, applying the operations of integration, wavelet 
decomposition, threshold operation, reconstitution and differentiation to the original signal 
5 by using the processing means of the computer, and deriving the main signal. 

Suppose a(t) is the original signal having a pseudo signal. The original signal is an 
acceleration signal detected by piezoelectric- type acceleration sensor 13 in the impact test 
as shown in FIG. 2A. First, original signal a(t) is integrated to obtain velocity signal v(t) 
(stepSlOl). 

10 Next, wavelet decomposition is applied to velocity signal v(t) to obtain the wavelet 
coefficient. The parameter for the wavelet decomposition can be, for example, Daubechies 
N = 3 for mother wavelet. In this connection, the procedure in which removal of unwanted 
component from the wavelet coefficient is followed by reconstitution to obtain the 
constituting signal of the original signal is hereinafter called "wavelet operation." Velocity 

15 signal v(t) is decomposed into a pseudo signal component (what is called Approximation in 
wavelet transformation) A(t) and main signal component (what is called Detail in wavelet 
transformation) D(t) by a predetermined wavelet decomposition (step SI 03). Then, a 
predetermined wavelet operation is carried out. Here, the predetermined wavelet operation 
is such an operation that satisfies following conditions. First of all, the energy ratio of A(t) 

20 to original signal v(t), E = || A(t) || V \\ v(t) || 2 , is calculated. Then, the A(t) which makes 
the energy ratio E smaller than a certain threshold value is derived, and the threshold level 
number which is the minimum of the level numbers as such is derived. In this way, by the 
wavelet decomposition at such threshold level number, D(t) of as many as the threshold 
number and one A(t) are obtained. In this regard, this threshold value represents the 

25 threshold value for the energy ratio of the pseudo signal component to the original signal, 
and it can be, for example, 99.7%. This is because, according to the method of the present 
invention, in the velocity signal derived by integrating the acceleration signal as the 
original signal, the pseudo signal component dominates large energy since the main signal 
does not have a zero-shift component and thus becomes almost zero when it is integrated. 

30 As a result, such A(t) can be regarded as a pseudo signal component. In this connection, 
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the predetermined wavelet operation just described is an approximation of that A(t) at the 
threshold level number to zero by the threshold operation to be described afterward 
followed by a reconstitution of the signal by using A(t) and D(t). In this connection, it will 
do if the wavelet decomposition is carried out up to the level which is deep enough to 
5 derive the above-described threshold level number. 

The wavelet operation is carried out for these D(t) and A(t) (step SI 05). The wavelet 
operation consists of removal of the pseudo signal component and reconstitution of a main 
signal. For the removal of the pseudo signal, first of all, there is a method where A(t) in the 
wavelet decomposition of the threshold level number is removed (made zero) as a pseudo 

10 signal component. In this case, residual D(t) represents the main signal. The wavelet 
reconstitution is carried out for the signal from which the pseudo signal component has 
been removed, and the resulting signal is output as velocity signal v c (t) corresponding to the 
main signal component. When A(t) is removed as a pseudo signal component, v c (t) is 
obtained by reconstituting D(t). On the other hand, the pseudo signal component is 

15 reconstituted by wavelet and is output as v^t) corresponding to the pseudo signal. When 
A(t) is removed as a pseudo signal component, v^t) is obtained by reconstituting A(t). 

Acceleration signal a c (t) corrected (from which the pseudo signal component has been 
removed) is obtained by differentiating v c (t), which corresponds to the main signal 
component obtained by removing a pseudo signal component A(t) from the original 

20 velocity signal v(t), with time (step SI 07). Thus, acceleration signal a c (t) from which the 
pseudo signal consisting of a zero-shift component has been removed is obtained. In this 
connection, pseudo signal a err (t) is obtained when velocity signal v crT (t) represented by the 
pseudo signal is differentiated with time (step SI 09). 

According to the characteristics of wavelet transformation, for the mother wavelet and 

25 its level, other mother wavelet or its level can be chosen depending on the characteristics of 
the target signals. 

[Working Examples] 

As shown in FIG. 2A, impact response acceleration signal a(t) having a pseudo signal 
measured by piezoelectric-type acceleration sensor 13 was used as an object of the method 
30 for processing a signal of the present invention. In this connection, the method for 
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measuring dynamic response characteristics by an impact acceleration meter disclosed in 
Japanese Patent publication for opposition purpose (Kokoku) No. 6-52270, in which a 
pseudo signal due to the zero-shift does not arise theoretically, was used as a comparative 
example, and the impact response acceleration signal was measured by a strain gauge as 
5 shown in FIG. 2B, thus the result obtained was used as a standard for quantification of the 
effect of the pseudo signal removal according to the present invention. 

FIG. 3 is a diagram showing impact response acceleration signal a(t) from which a 
pseudo signal is to be removed. Velocity signal v(t) (FIG. 3) was obtained by integrating 
impact response acceleration signal a(t) shown in FIG. 3. Each of a(t) and v(t) has such 

10 shape shown in the figure because of the effect from the reflected wave at the end surface 
of the rod. Then v(t) was decomposed into main signal component D(t) (FIG. 5) and 
pseudo signal component A(t) (FIG. 6) by using wavelet decomposition. In this example, 
Daubechies N=3 was used as the wavelet. In this connection, threshold operation was 
carried out by assuming the energy ratio E of pseudo signal A(t) to the original signal v(t) 

15 was 99.57%. For this purpose, the wavelet transformation is executed with deepening the 
level number until the threshold level number at which the ratio of the energy E x has 
predetermined relationship with a predetermined threshold value. Namely, A(t) obtained 
by the wavelet decomposition at the threshold level number which was the minimum of the 
level numbers where the ratio of the energy E became smaller than 99.51% was removed as 

20 a pseudo signal component. Main signal component D(t), from which pseudo signal 
component A(t) has been removed, and which corresponds to corrected velocity signal v c (t), 
was differentiated with time to obtain corrected acceleration signal a c (t) (FIG. 7). In this 
connection, acceleration signal of the pseudo signal component a eiT (t) (FIG. 8) was obtained 
by differentiating pseudo signal component A(t) with time. FIG. 9 is a diagram showing 

25 the results of Shock Response Spectrum (SRS) analysis for a(t), a c (t), a crr (t) and impact 
response acceleration signal a^t), which was deemed as true value, measured by the strain 
gauge according to the comparative example. In FIG. 9, each SRS was denoted by "O" for 
a(t), "+" for a c (t), for a^t) and "A" for a CIT (t), respectively. Referring to FIG.9, it is 
observed that a c (t) fits a^t) well, which means that the pseudo signal component was 

30 removed by the method of the present invention. 
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From the following facts, it is understood that the precision of the zero-shift removal 
was improved by the method of the present invention, in comparison with the prior art. 
FIG. 10 is diagrams showing the original acceleration signal having a zero-shift (upper 
diagram) and its SRS (lower diagram). Referring to the acceleration signal, although it 
5 converges from the initial oscillation having large amplitude to a constant value, the center 
thereof deviates significantly upward from zero at any point of the time period, which 
shows that it has a zero-shift. In addition, referring to SRS, it is observed that the positive 
component denoted by "O" significantly deviates from the negative component denoted by 
" A ", and is clearly recognized that the acceleration signal deviates to the positive side. 

10 FIG. 1 1 is a diagram showing the acceleration signal (upper diagram) which is an analytical 
result obtained by carrying out the wavelet operation directly (without carrying out the 
integration) to the acceleration signal according to the prior art method and its SRS (lower 
diagram). Referring to the acceleration signal, although it converges from the initial 
oscillation having large amplitude to a constant value, the center thereof deviates slightly 

15 upward from zero especially at the beginning, which shows that it has a zero-shift. In 
addition, referring to SRS, it is observed that the positive component denoted by "O" does 
not fit the negative component denoted by " A and is also clearly recognized that the 
acceleration signal deviates slightly to the positive side. FIG. 12 shows the results obtained 
by integration of the acceleration signal to transform it to the velocity signal and to 

20 emphasize the zero-shift component followed by wavelet operation according to the 
method of the present invention. Referring to the acceleration signal, it converges from the 
initial oscillation having large amplitude to a constant value, and the center thereof locates 
near zero at any point of the time period, which shows that a zero-shift is scarcely observed. 
In addition, referring to SRS, it is observed that the positive component denoted by "O" 

25 nearly fits the negative component denoted by " A " throughout the whole region, and is 
clearly recognized that the acceleration signal has no deviation due to the zero-shift, which 
means that the pseudo signal component has been removed effectively. From these facts, it 
is understood that the precision of the zero-shift removal was significantly improved by the 
method of the present invention, in comparison with the prior art. 



13 



As described above, according to the present invention, the pseudo signal component 
can be effectively removed from the original signal which has the pseudo signal component 
consisting of zero-shift component. 
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